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CMG-INDUCED LST DYNAMICS 
INTRODUCTION 


The nature of the control moment gyroscope (CMG) output gimbal bear- 
ing friction has been the subject of many discussions and presentations by 
various investigators. It is important to know the precise nature of the torque 
caused by the bearing friction so that its effect on the dynamic behavior of the 
Large Space Telescope (LST) and its pointing and control system (PCS) can be 
predicted. The objective of the analysis and supporting computer simulations 
reported herein is to study the effect of the assumed CMG nonlinearity upon 
the system dynamics in order to develop an understanding of the nature and 
effect of variations of numerical values of parameters upon the system. The 
weakness of the mathematical analysis used results from the need to restrict 
the complexity of the model studied; however, the design insight obtained is 
valuable. 

The parameter plane analysis technique used herein makes use of the 
describing function to portray the nonlinearity representing CMG gimbal bearing 
friction. Where the describing function is a complex value, the technique 
developed for a system containing two nonlinearities may be applied [ l] . This 
technique affords both analytic and graphic portrayal of the effects of variations 
in selected system parameters. With few exceptions (noted herein) the results 
are confirmed with analogue computer simulation. 

Finally, the reader is cautioned against relying solely on computer 
simulations without adequate mathematical analysis as a background. An 
attempt to use an analogue computer to define stability boundaries for satellite 
dynamics as characterized by the Mathieu equation should convince even the 
skeptic of the weakness of computer simulation when not augmented by mathe- 
matical analysis. A current example in the field of digital simulation is the 
errors contained in Connell's paper (2) reported by Schiehlen [3]. 

The following four models are considered in this report; three of them 
are analyzed in detail. 


1. A model developed by the Bendix Corporation and presented to MSFC 
in April 1972 [4-6 ] 1,2 . A single axis representation of the CMG, rate gyro, 
and rigid body dynamics is shown in Figure 1. 


CMG 



Figure 1. Bendix model number 1. 


2. A Sperry model (Fig. 2) presented to MSFC in December 1972. 

3. Another Bendix model presented in March 1973 (Fig. 3). 


4. The model developed by Dr. P. R. Dahl of Aerospace Corporation 
( the so-called Dahl model) . This model may be incorporated in a single axis 
representation (Fig. 4), where the nonlinear relation between the frictional 
torque (T^) about the CMG output gimbal pivot and the gimbal rotation (6) 


is portrayed as the nonlinearity N. 
by Dahl, of this relation is 


T 


f 


y5 


» 


The mathematical description, develooed 






(i) 


1. CMG Considerations for HEAO and High Accuracy Point Missions. Written 

Presentation Material, Bendix Corp. Teterboro, N.J. , April 21, 1972. 

2. Whitley, G.W. : LST Fine Pointing Control System Design. MSFC S&E- 
ASTR-SD-78-72 Letter, Oct. 3, 1972. 
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and is shown in Figure 5, where y and T GFq are parameters describing the 

friction relationship and the overdot represents differentiation with respect to 
time. 



Figure 5. CMC friction nonlinearity a la Dahl. * 

m 

ANALYSIS TECHNIQUE 


The objective of the analyses presented herein is to determine stability 
''onditions tor systems that may be represented schematically, as in Figure 6. 
The niant under consideration is represented by a linear part G(s) and a non- 
linear elei .ent N. In the sequel N will ne used to represent the nonlinear 
friction characteristics of the CMG gimbal pivot. In particular it is of interest 
to Investigate the possibility of existence of limit cycles and to determine their 
characteristics. As used herein the term "limit cycle" will be used synono- 
mouslv with the term "self-excited sustained oscillation." Hence, a limit 
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Figure 6. Basic r. nlinear system. 

cycle represents. . steady -state oscillation to which or from which all neigh- 
boring trajectories in a state space will converge or diverge [7,8] . The 
stability of limit cycles may be defined in terms of three classes: 

1. Stable ( Orbitally Stable) Limit Cycle — If all neighboring state 
space trajectories converge to the limit cycle. 

2. Unstable (Orbitally Unstable) Limit Cycle — If all neighboring 
s^ate space trajectories diverge irom the limit cycle. 

3. Semistable (Orbitally Semistable) Limit Cycle — If state-space 
trajectories originating at points outside the limit cycle converge to it and 
thop~ originating from points inside the limit cycle diverge from it, or vice 
versa. 


Because of the low pass nature of the linear parts of the system plants 
under consideration, they are particularly amenable to analysis through 
application of describing function theory. Hence a detailed discussion of 
aporopriate describing functions is needed [9l . It is assumed that, if a 
sustained oscillation occurs, the input x is described by 


x = A sin ftt . (3) 

T t is assumed that the output f(x) is approximated by the relation 

f(x) » NA sin (ft + <f>)t , (4) 


where $ represents a phase shift. 
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The method o£ analysis used is defined in Reference 1 and embodies a 
describing function technique developed by Siljak [ 10] . It is assumed that the 
system being analyzed is amenable to describing function analysis and the 
system parameters are time-invariant. In essence, a describing function (N) 
is used to represent the nonlinearity. The characteristic equation representing 
the system is then obtained. For a limit cycle to occur, it is required that all 
characteristic equation roots have negative real parts except for one pair of 
roots which must be purely imaginary roots. This condition is determined 
mathematically and for ease of visualization may be portrayed graphically. In 
the latter case, two adjustable parameters are selected, at least one of which 
must also contain N. (if N is a complex quantity, its real and imaginary 
components may be used as the two parameters and Reference 1 may be applied 
directly.) A correlation between these parameters and the roots of the charac- 
teristic equation is then determined by mapping stability contours from the 
complex s-plane onto the chosen parameter plane. For the systems considered 
herein, only three (and sometimes fewer) such stability boundaries exist and 
they are easily found. 

The first boundary (the one associated with real roots of the charac- 
teristic equation) is found by setting s = 0 in the characteristic equation. 
(Graphically this boundary may be identified by single hachures. ) The second 
one (the one associated with purely imaginary roots) is found by setting 
s = icu . (This boundary is identified by cross-hachures. ) The third one, called 
the boundary at infinity (in topological mapping from the complex s-plane to the 
Riemann sphere, it is associated with the sphere's north pole), is found by 
setting l/s = 0 in the characteristic equation. It also is identified by single 
hachures. From these boundaries, the stable region (if it exists) may be 
determined in terms of the selected adjustable parameters, for they bound the 
stable region. 

The nonlinear locus of N as a function of the two adjustable parameters 
is determined next. The simultaneous solution of the nonlinear locus relation 
and the purely imaginary root boundary yields the condition for a limit cycle, 
assuming the indicated solution occurs adjacent to a stable region (which may 
not always be the case). This is readily apparent on the parameter plane as 
the intersection between the purely imaginary root boundary and the nonlinear 
locus. From this point of intersection, the frequency (identified as = fi) and 
magnitude (A) of an indicated limit cycle may be determined as a function of 
the characteristics of the nonlinearity and of the adjustable parameters. 

Further, the behavior of the limit cycle when a small perturbation is applied to 
its amplitude, and, nence, the nature of limit cycle stability, also is apparent 
on the parameter plane. 
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ANALYSIS OF BENDIX MODEL NUMBER 1 (Fig. 1) 


Four variations of Bendix model number 1 were studied and are identified 
herein as Cases I through IV: 

1. Case I — This model consists of a perfect attitude rate sensor 
characterized by 


R(s) = 1 


(5) 


and a CMG characterized by the nonlinear element N and a first order trans- 
fer function 


°g< 8 > = 


CD 

c 


S + CD 

c 


( 6 ) 


This is the basic model considered. 

2. Case II — A model characterized by a perfect attitude rate sensor, 
equation ( 5) , and a CMG modeled as a nonlinear element N and a second 
order transfer function 


G g (s) 


K. 


+ gj S + g Q 


(7) 


The effect of a higher order model of the CMG is induced by comparing the 
results of the analyses of Cases I and II. It is similar to several models con- 
tained in the documents cited in References 4 and 5 and footnote 2. 

3. Case III — A model consisting of an imperfect attitude rate sensor 
described by a first order transfer function 


R(s) 


CD 

£_ 

S 4 CD 

K 


( 8 ) 


8 




r 


1 


and a first order characterization of the linear portion of the CMG, equation 
(6) . This model permits an analysis of the effect of a nonperfect attitude rate 
sensor. It is the model documented in Reference 6 and footnote 2. 

4. Case IV — An imperfect attitude rate sensor, equation (8), and a 
second order characterization of the linear portion of the CMG, equation (7) . 
This model is similar to several models documented in References 4 and 5. 

Because the precise nature of the relationship between the input (x) to 
the nonlinear element (N) of the CMG model and the output f(x) is obscure, 
the effect of variations of the parameters describing that relationship is studied. 
In particular it is desired to predict whether or not sustained oscillations in 
the vehicle's attitude and attitude rate will occur and, if so, what their charac- 
teristics are ( in terms of amplitude, period, and stability) under various 
prescribed conditions. 


Basic Nonlinearity 

If the small overshoots occurring at x = D and x = - D (Fig. l) are 
ignored, the describing function is 


7 " = 1 


- 1 [(* - £) sin '' u + (* - f - t) "' JTr:? ] 


, u < 1 
(9a) 


= m , u > 1 


(9b) 


For small values of m (m« k), the relation approximated by the describing 
function defined in equations (9) then becomes 




- £sin -1 u + ^1 - 1 -u 5 j , u ^ 


(10a) 


= 0 , u > 1 


(10b) 


9 


i 


1 


whore 


u 


d. 


D 

A 


(10c) 


and 


d. Dk 
M 


( iod) 


Although the describing function may be found in a number of texts on nonlinear 
ose illations, this author was unable to locate a sketch of the variation of the 
describing function N with variations in the parameters A , 1), M, and k. B\- 
selccting combinations of these parameters in a certain way, an interesting 
analogy between this particular describing function and a system familiar to 
control system engineers was discovered. Symbols i), M, and k are described 
on Figure 1; I) represents the dead zone, and M and slope k characterize the 
gam through N. If N k is plotted versus A D with £ as an independent param- 
eter (Fig. 7), it is observed that all curves begin from a common point at 

N 0 tor A l) i and asymptotically approach a value of unity for large 
values of A I). If (2 1 the curves never rise above the value of unity. How- 
ever, if is 1 the curves rise from their initial value of zero to a peak value 
greater than unity and then (as A/D increases) approach the unity asymptote 
tiom above. The analogy to a second order control system is interesting. 

Here, 5 is analogous to the damping ratio 4. Further study of the describing 
function reveals the amplitude of the peak value of N ''k and the values of £ and 
u (and hence D, M, and k) for which the peak occurs. If equation (10a) is 
optimized with respect to u, it is found that the peak value, N n k, of N k 
occurs when u has the value 



00 


t.eneral (i.e., independent of numerical values for A, I), and k) curves of 
N o /k versus < and A ft /D ( l/u n ) versus ( may be plotted from equations 


( 10a) and ( 1 1) and arc shown in Figure S. These relationships will be used 
in subsequent predictions of limit cycle conditions. 
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Figure 7. Describing functions for Bendix model number 1 
nonlinearity ( simplified) . 


Additional insight into the nature and dynamic effect of the nonlinearity's 
parameters A, D, and k may be obtained by considering certain limit cases. 

Limit Case 1 (Ideal Relay) . If the dead zone (D) of the basic non- 
linear ity(FTgTTy7snegiected^nd slope k decreased to zero, an ideal relay 
characteristic is obtained where 


N = 


4M 

-A 


( 12 ) 




A diagram of this relationship is shown in Figure 9, setting k = 0 . 

Limit Case 2 ( Relay With Dead Zone) . If the dead zone of the basic 
nonlinearity is retained but slope k set at zero, a relay with dead zone charac- 
teristic is obtained where 


N 


4M \' 1 - u- 
-A 


u ^ 1 


(Li) 


This relationship is shown in Figure 10. Observe that DN/M has a peak value. 
This may be found by setting $ = 0 (since k - rt) in equation ( 11) and sub- 
stituting Up = D/A = lA l~2 in equation (12) yielding a peak value of DN/M = 

2/«r « 


A. 


J 
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Figure 11. Describing function relation for gain with dead zone. 

Limit Case 4 (Preload) . If the dead zone of the basic nonlinearity is 
collapsed to zero, a preload characteristic results, yielding 


H. « ii + i 

k 


This relationship is shown in Figure 9 with k # 0 . 

Limit Case 5 (Gain Change) . If the parameter of M of the general non- 
linear ityTFig^n~^s~setTquaTlo^ero, a gain change characteristic results, 
yielding 

j“ = 1 - ~ y) ( sin_1 u + 1 - u 2 ) , u < 1 , (16) 

= 1 , U s 1 . ► . 

The relation is shown in Figure 12. 

The four cases will now be analyzed in detail, considering the above 
limit cases [equations (12) through (16)1 and the general case modified by 
setting m equal to zero [ equations ( 10) 1 . 


Analysis of Case I 


The characteristic equation is 


s 3 + co s 2 + BNs + AN = 0 , 

c 


14 



Figure 12. Describing function relation for gain change. 


where 


a £• JLs 


j K« 

B =• -i_£ 

V 


It is readily shown that the stable region of operation, in terms of 
parameters^ AN , BN lies in the first quadrant and is bounded by the real root 
boundary (AN =0) and the imaginary root contour 


BN = 


From definitions ( 18) , the nonlinear locus is seen to be a straight line 
through the origin with slope K 1 /K q . These two lines are shown on Figure 13. 

Since they cannot intersect, limit cycle operation is not indicated. A stable 
response ie indicated if K /K > l/u> and, if K /K < l/u> , the system 

A U C 10 C 

response will be unstable. 


15 


t 



Figure 13. Case I (Bendix model number l) — stability characteristics. 

If the pitch axis is chosen as the single -axis about which rotational 
dynamics will be studied, the numerical values given in Table 1 may be used 
as representatives of the LST. With the exception of co and , these values 

c g 

were obtained from Reference 5. The value for is obtained by matching 

the break frequency o> c [equation (6)1 with the break frequency of the second 

order system: 


W c ' v ' g(/ 2 ~ 70.7 rad/s 


( 20 ) 


The value selected for is representative of a bandwidth-limited rate 

gyroscope and is low enough to indicate the dynamic effect of this degradation. 

Using Table 1 values, it is seen that a stable response is indicated, i.e. , 

K /K > 1/co . 

10 c 


Analysis of Case 1 1 

The characteristic equation is 

s 4 + g^ 3 4 g^s 2 s ^ Ns 4 k N = 0 (21) 


♦ 


t 


16 




Parameter 

Numerical Value 

g 0 

= K = 10 4 (rad/s) 2 
G 

g l 

141.4 rad/s 

H 

610 Nms 

*GE 

5.012 Nms 2 (3.7 ft lb s 2 ) 

I 

V 

1.356 x 10 s Nms 2 

K o 

4.6843 x 10 6 Nm 

K 1 

1.1153 x 10 G Nms 

k F 

6.600 x 10 3 Nm/rad 

^CMG 

0.707 

K 

0 

4.352 x 10 4 (= oj 2 I - b.) 

n G t 

K 

1 

7.087 x 10 4 (= 2 ; o) I ) 
n G 

C 0 

c 

70.7 rad/s 

Cl) 

K 

30 rad/s 

W ] 

n CMG 

100 rad/s 








1 


where 


k =’ K K /I , 

0 0 G v ’ 


(22a) 


S - K , Vv 


(22b) 


The real root boundary is 


k N = 0 
0 


and the imaginary root boundary, in terms of parameters k N and k IS , may 


be written as 


k o N - 


/ k l N 

k i N ( *o - ^ 


or, written as a function of the independent argument, frequency (cu). 


k Q N = cu 2 ( g() - cu 2 ) 


(25a) 


k,N = g t u) 2 


(25b) 


From equations (22) , the N-locus is expressed as 


S N ’ E 


A limit cycle and its characteristics are Indicated by the simultaneous solution 
of equations (25) and (26), yielding 


Vv ( 

N “ k i k g V 0 


K o g i 


Si = v g. 


Vi 


The trequency n of the Indicated limit cycle Is determined by system pa ram- 

A M J . 11 euf firrua rtf nrtn _ 


I HC Li ^^UVlIVJf ““ V* vaa-*" V 

eters g • g . K • ^ and has the same value regardless of the type of non- 
linearity chosen for N . Although the value of N is also independent of the 
type of nonlinearity, the amplitude A of the limit cycle is dependent on the 
relation between N and A, e.g., equations ( 10) , and (12) through (16). 
Equations (23) , (25) , and (26) , are plotted on Figure 14. The simultaneous 
solution of equations (25) and (26) is indicated by the intersection labeled 
limit cycle. 

Before numerical values are used, some general characteristics may 
be deduced. If the nature of the nonlinearity is such that N increases rrono- 
tonically as A increases (curve ( 1) , Fig. 14), then an unstable limit cycle 
is indicated. However, if N increases but approaches a limit asymptotically 
as A increases, the limit cycle wiH occur only if the asymptote has a value of 
k N corresponding to the limit cycle (for curve (2), the value is less, so 

no limit cycle is indicated] . If N decreases monotonically as A increases, 


r 



stability characteristics 



a stable limit cycle will occur as long as the curve of k^N versus A begins 

at a value greater than k N [e.g. , curve (3) ] . If the curve remains above 

1 LL 

the value of k N for all A [curve (4) ], no limit cycle will occur. These 

1 JJV 1 

observations lead to an investigation of describing function characteristics for 
various appropriate nonlinearities. 

Limit Case 1 (ideal Relay) . Since N decreases monotomcally as A 
increases, a stable limit cycle is indicated with an amplitude [see equation (12) 1 

A = 4M/Nir . (29) 

[N is obtained from equation (27) and frequency 12 is obtained from equation 
(28). I 

If the numerical values of Table 1 are used, 


10.1707 

(30) 

90.985 rad/s = 15.4357 Hz 

(31) 


§ 

If, for example, M is chosen as 0.055, a value of A = 0.00434 is 

predicted. 3 The amplitude NA of the assumed sine wave f(x) is NA = 

0.07018. The amplitude of I 0 1 and 1 0 1 of oscillations in 0 and 0 

LC LC 

may be determined from this relationship with the signal x (Fig. l). Assum- 
ing the form of equation (3) for x, the amplitude of the 6 limit cycle 

and the 0 limit cycle 1 0 1 may be approximated by 

LC 


I 0 1 


LC 





(32a) 


3. Superscripts u and s refer to quantities associated with unstable and 
stable limit cycles, respectively. 


21 


and 


101 


LG 



(32b) 


I’sing equations (32) , amplitudes of 1 0 ! and Ifll are found to be 

La. LiL 

4.0123 n 1 0 ~ ’ rad and 3.891 x 10~ 9 rad's, respectively. 

Limit Case 2(Relay With Dead Zone) . In this ease stable operation 
without a limit cycle is indicated if DN M v 2 / 'tt . If |)N M' 2 ", both a 
stable limit cycle ami an unstable limit cycle are predicted, since the curve of 
kjN versus A rises above the value of k ^ N ^ to a peak value and then falls 

below it. A case of academic interest arises when 1)N M is precisely equal 
to 2 " . In that case, a semistable limit cycle will occur. F sing a value of 
M 0.55 and l) 0.0011 , equation (13) may bo solved for the two predicted 

values of A, A U 0.00114 and A‘ S 0.0042 , and amplitudes of f(x), NA U 
s 

0.0184 ami NA^ 0.0(579 . From equation (32) , 0 and 0 limit cycle magni- 
tudes are four * to be |0|^ 1.0539 x 10~" rad, |0| ^ 1.0221 ^ 10 -9 rad/s 

1A. L(v 

Limit Case 3 (Gain With Dead Zone) . It is seen from equation ( 14) 
and Figures 1’ and 14 that an unstable limit cycle will occur if k -• N , and 
stable operation will occur if k N . An example is obtained by setting k 20 
If I) 0.011 , equation ( 14) may be solved for the predicted value of As 

A U 0.0729. From equation (32) , I <» I V (5. 73% ^ 1 0* to rad and Id 1 11 

lA l A 

(5. 53(54 N 10" rad s. Amplitude of f(x) then becomes NA U 1. 1788. If I) 
0.00022. then .\ U 0.0014(5, NA 0. 023(51, Id I* 1 1.3498 x lo -M rad, and 

I (> I j 1 . 309 1 x 1 0' 9 rad s. 

Limit fast' 4 (Preload) . It is apparent from equation (15) and Figures 
9 and 14 that a stable limit cycle is predicted if k v N and that unstable opera- 
tion will occur if k ^ N . An example of the effect of a preload characteristic 
nonlinearity is obtained by setting k l , M 0.055. Solving equation (15) fot 

A, one predicts A S 0.004(5 111 ami N.\ S 0.074(54 . From equation (32) , 

Itflj 4.2(57 v to-" rail and Idl* 4. 1388 \ 10" 9 rad s. 


the value of asymptote N/k~ 1 when 4 > 1 ( k > M/D). Using the same means 
of analysis as used in the preceding cases, numerical examples corresponding 
to each of these six conditions are developed. They are summarized with the 
results of the numerical examples of the five limit cases in Table 2. Attention 
is directed to the case covered by Condition 1 where M = D , k = 1 : This is 
the set of conditions on ‘he nonlinearity that is used in References 4 through 6 
and footnote 2 and is reported as causing limit cycles. The analysis contained 
herein disagrees with the reported results in that this analysis predicts stable 
operation and the absence of any limit cycles. 


Analysis of Case III 

The characteristic equation is 


S 4 + (u) S 3 

C g 


c jv 


K u) 4 k ) - — 
C l £ 0 1 

H V 


I 


^0 

R 



v 


0 


9 


(34) 


where overbars are used on gains and to identify them with the system 

of equation (34) rather than Case II. Comparing equation (34) to equation (21), 
one can immediately draw conclusions from the analysis associated with the 
system represented by equation (21) . Matching the terms. 


g = u) •+ a.' 

leg 


g . r U? 

0 c g 


k - K /I 
0 0 0 v 


and 


k, u) + K ) /I 
1 c 0 g 1 v 


(35a) 

(33b) 

(35c) 


( 35d) 
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Limit Case 5 (Gain Change ). For the anticipated values of m « k , 
equation (Hi) is approximated by equation ( 14) . It is anticipated that the 

results ot the analysis of Limit Case 3 will apply approximately to Limit Case 
5. 


Basie Nonlinearity (Fig, l) . The effect of varying the basic non- 
linearity s parameters k, M, and D will be investigated in detail. An examina- 
tion of Figures 7 and 14 yields most of the following information: 

Condition l. N > k r* M 'D: Stable operation (no limit cycle). 

Condition 2. N > k < M D - k'^s Two limit cycles; one stable, one 
unstable with A u A‘ s . 

( ondition :t. N'kv M D v k Stable operation ( no limit cycle) . 

Condition 4. N k - M D: Unstable limit cycle. 

Condition 5. N k v M l): Unstable limit cycle. 

Condition (>. N ■> k v M I) k 'i. : Semistable limit cycle. 

Additional information was needed to obtain the inequality relationships 
between M l) and k for conditions 2 , 3 , and ti. From equation ( 11) , it 

is seen that a peak in N k versus A I) can occur only when 


£ 


(33a) 


or 


( :i:tb) 


Since l) A must be less than unity for equation (10a) to bo applicable 
it is seen from equation ( ll) that inequality Ckta) is meaningless physically. * 
Also from equation ( 10b), it is seen that • - 0. To summarize, N K will 
have a peak only when IK > , 1 (k . M/D) and will approach but not exceed 
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SUMMARY OF BENDLY MODEL NUMBER 1 RESULTS 







































































































































Now, equations (21) through (26) may be used directly for this model. 
Hence, it is seen that the effect of — i. e. , a nonideal rate gyro — is to 

introduce the possibility of a limit cycle, considering Conditions 1-6 listed 
previously. A numerical example is selected using the values in Table 1, 
yielding N=9.69 and fi = 40.64 rad/s. If a nonlinearity of M=D and 
k = 1 is chosen, it is again seen (Condition l) that limit cycle operation will 
not occur. The analysis used in Case II and summarized in Table 2 is directly 
applicable to this case as well. 


Analysis of Case IV 

The characteristic equation is 


s 5 + (gj + u^)s 4 


t m aN = 0 , 

g 


(s. 




CO g s‘ 

go 


+ ( u j /3 + cv) Ns 

ft 

S' 


(36) 
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r 


where 


a 


d. 


K K /I 

0 G v 


(37a) 


and 



K , K /I 
1 G v 


(37b) 


The stable region, in terms of an evN -d N parameter plane, is bounded by 
the real root boundary 


«N = 0 


(33) 


» I 


♦ 
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and the imaginary root boundary 


I 





(39a) 


and 


0N = g co 2 (l + u> 2 /a ) 2 ) 
1 & 


(39b) 


Again, the nonlinear locus is the straight line determined by equation 

(39): 


0N = (K /K ) 


(40) 


This is shown in Figure 15. If a limit cycle is to occur, it will have the 
characteristics associated with the intersection of the boundary defined by 
equations (39) and the nonlinear locus [equation ( 40) ] - Using the numerical 
values of Table 1 to solve equations (39) and (40) simultaneously leads to 


N = 2.103 


(41) 


and 


ft = 40.959 rad/s = 6.519 Hz 


(42) 


The six conditions listed previously may be used to analyze the possi- 
bility of limit cycle operation. Again, it may be readily concluded that for 
M = D and k = 1 , limit cycle operation will not occur. Limit cycle charac- 
teristics for other values of M, D, and k may be predicted, however. Examples 
are summarized in Table 2. 
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irst quadrant view 












Most of the results presented above for Bendix model number 1 have 
been confirmed by analogue simulation. 4 The model of Figure 1 was set up on 
the EAI 231R-V computer in the Astrionics Laboratory at Marshall Space Flight 
Center. The simulation diagram used is shown in Figure 16, and the nonline- 
arity used on the computer is shown in Figure 17. This simulation description 
is detailed and lengthy because it graphically portrays numerous combinations 
of unusual limit cycles normally found only in textbooks. 



Simulation of Case I 

Several simulation runs were made to confirm the predicted stable 
behavior of the system. The simplicity of the results does noi warrant repro- 
duction herein, so they are omitted. 


Simulation of Case II 

This case provides the basis of the results reported herein an^is 
described in detail. The tabulated results (Table 2) will assist the reader in 
following the discussion. 

4. The results presented in this section are due to the efforts of Mr. P. H. 
Fisher, Guidance and Control Systems Branch, Systems Division, Astrionics 
Laboratory, MSFC, who worked with the author to confirm the analytically 
derived results. His contributions are gratefully acknowledged. 
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Figure 17. Analogue simulation diagram for CMG nonlinearity of 
Bend Lx model number 1. 

Limit Case 1 (Ideal Relay) . A stable limit cycle is predicted and is 
confirmed in Figures 18 and 19. In Figure 18, the limit cycle trajectory is 

approached from within (in a state space sense) with x(0) < A" , and in Figure 

g 

19 the limit cycle trajectory is approached from the outside with x(o) > A . 
Hence, the stability of the limit cycle is demonstrated. The steady state 

amplitudes of x, i.e. , A S , 0 and 0 , may be read directly from Figures 18 
and 19 and are tabulated in Table 2 for comparison with the predicted values. 
Only the value of the magnitude of f(x) varies appreciably from predicted 
values because of its assumed (for describing function analysis purposes) 
sinusoidal shape. The frequency of the sustained oscillation closely matches 
the predicted frequency ( fi) in this case and all other Case II simulations. 

Limit Case 2 (Relay With Dead Zone). Both the unstable and stable 
limit cycles are predicted for this case. The amplitude of A was too small 
(i.e. , too nearly identical to dead zone D) to be observed. However, the 
stable limit cycle was detected and its nature confirmed. A sample run, with 

xi •))m' A^, is shown in Figure 20. 

Limit Case 8 (Gain With Dead Zone) . Because of the close agreement 
between analytical and simulation results in other cases, this case was not 
simulated. 


30 





0 . 5*10 9 r 


0.5* 10- 9 


0 . 5*10 9 



0 1 OJ 


Ki^urt' 19, Limit Case 1 (Ideal Relay) simulation, M 0.095: 

\( 0 ) ' A S . 

Limit Case 1 ( Preload) . This is another case of a predicted staid 
limit cycle. Similar to the approach of Limit Case l, the predicted staldi 
limit cycle characteristics are confirmed in Figures Jl anil 22 . 
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I' *£ uro I-imit Case 4 (Pro loud) simulation, k l, M 0 . 0 ;>;>: 

x(0) - A 8 . 

Limit t aso f> (« ’. ain i-hango) . This case, similar to Limit Case :l since 
m w k, was not simulated. 
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Condition 2{N>k<M/D> k/4 Q ) . In this condition (as in Limit Case 

2) both unstable and stable limit cycles are predicted. As in the nature of 
unstable limit cycles, it is extremely difficult to maintain their oscillations. 
However, Figure 24 contains a good example of an unstable limit cycle main- 
tained for several periods before the trajectory escapes to the sustained 

oscillation of the stable limit cycle (A > A ). In another run (not included in 
report) , the unstable nature of the limit cycle was seen as the states decreased 
to zero. The stable nature of the larger amplitude cycle is confirmed in Figure 

25 where x(o) > A b . 


t 

e 



-OJnlO-^h 



Figure 24. Basie nonlinearity simulation, k- 4.95, 13 0.011, 

M - 0.275: x (0) < A U . 
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Figure 25. Basic nonlinearity simulation, k - 4.95, I)- 0.011, 

M - 0.275: x (0) - A 8 . 

Condition :i (N k ■- Md) v ). As predicted, stable operation 
ensues ( Fig. 20 ) . 





Figure 26. Basic nonlinearity simulation, k= 0.95, D= M 



Condition 4 (N - k > M/I)). This ease (Fig. 27) provides an unusually 
iv ( )od example of an unstable limit cycle maintaining a sustained oscillation lot 
a large number of periods before (in this case) the response diverges from 
the limit cycle trajectory. 



0.1 0-2 


Figure 27. Basic nonlinearity simulation, k 20, 
I) M 0.011 


Condition r> (N k . M D) . Figure 2S is an example of an unstable 
limit cycle. The first run shows the elusive unstable limit cycle trajectory 
rapidly diverging to a stable response. The second run shows the trajectory 
just as quickly, diverging to an unstable response. 






Condition 6 (N < k < M/D - K/£^). This run (Fig. 29/ confirms the 

predicted behavior of a stable limit cycle and an unstable limit cycle merged 

u s 

into one, yielding a semistable limit cycle. In this run, x(o) > A ’ and the 
system response approaches and maintains the limit cycle behavior. After a 
number of cycles elapse, however, the response leaves the limit cycle tra- 
jectory to assume a stable response. Hence, the limit cycle is stable when 
the trajectory is approached from outside and unstable when approached from 
inside. 
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Figure 29. Basic nonlinearity simulation, k= 1, D= 0.011, 

M = 0.257: x (0) > A U,S . 

Simulation of Case 1 1 1 

Because of the close relationship to Case II, this case was not simu- 


lated. 


Simulation of Case IV 


This case should have been confirmed by simulation but was concluded 
because of time constraints. 


Summary 


Because the describing function technique assumes that the input and 
output to the nonlinearity are both sinusoidal in nature, it is to be expected 
that the actual magnitudes of these transcendental functions will only be approx- 
imated by the mathematical estimates. The closeness of the predictions to the 
simulation results is surprisingly accurate, considering the shape of the wave- 
forms at the output of the nonlinearity. The predicted time periods were 
particularly accurate. 


ANALYSIS OF SPERRY MODEL IFig. 2) 


A brief (compared to the preceding analysis) analysis was made of the 
Sperry model that was presented to MSFC in December 1972. The character- 
istic equation associated with the model is 



P 
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where N represents the CMG nonlinear element (gimbal pivot friction). If 
l/k F and N are selected as the two parameters of interest in the limit cycle 

and stability investigation, one obtains three stability boundaries. The real 
root boundary is 


N = 0 


(44) 


and the boundary at infinity is 

* V 

l/k F = 0 


(45) 


The imaginary root boundary may be written as 


l/k F = & (I V I GE & . (I v k o + HK^) fi* + HK o k o 1/J (46a) 


N = \ “ 4 fI /i 02 - H(K o K l + K 1 K 0 ),/J 


(46b) 


where the -Jacobian J of the simultaneous equations obtained from the real and 
imaginary parts of equation (43) is 


j = - n 2 


[Vi * - »<Vi + Vo>] 


"v’o ' (I v K o + + HK 0 « 0 I : • 


( 46c) 


For the general case, the l/k p - N parameter plane stability plot 

takes the form of Figure 30. The region of interest of this figure for limit 
cycle considerations is enlarged and plotted on Figure 31. Shown in Figure 10 
is a sketch of the describing function relationship associated with the nonlinearity 
Is of the Sperry model (Fig. 2), where 


44 


N 



Figure 30. l/k - N parameter plane — Sperry model. 

F 

N = 4M\' 1 - (D/A) 2 /VA , A/D 2 1 . (47) 


If the nonlinear locus N (as a function of A) is reflected from the describing 

function graph onto the parameter plane, it is seen to be a vertical line rising 

from (and returning to) l/k as A increases positively. Thus if l/k > 

r r 

1. 16 x io ~ 7 (k < 8.62 x 10 6 ) , no limit cycle can occur. If l/k < 1. 16 x 10“ 7 , 
F r 

then /.ero, two, or four limit cycles are predicted, depending on whether the 

peak value of ND/M(= 2/ff) intersects with none, one, or both of the portions 

of the complex root boundary. From Table 1 it is seen that the nominal value 


45 


1 ^ 


N 



Figure 31. Limit cycle region of interest on parameter plane — 

Sperry model. 

of k is 6600 ( l/k., =-■ 1. 52 * 10“ 4 ) , indicating a wide margin of safety (for 
F r . 

the numerical value chosen for the system parameters of this model) tor 
precluding limit cycle operation. For example, if M - 10D- 2.2 ^ 10* and 
1/k , - h lO" 7 , then the lower value of N is 0.30 10~*' (with a corre- 

sponding value of 0 -■ 0.57) and the upper value is 1.95 x 10" 6 (and SI « 1.06) . 
Since l) is so small (with respect to unity), the two values of A corresponding 
to each value of N are [from equation (47)1 


A « 4M/*N , 0 


(4S) 


>6 


Thus two of the four indicated limit cycles have negligible amplitudes (A) and 
can be ignored. The other two significant amplitudes, the one corresponding 

g 

to N = 0.3 x 10 -6 being a stable limit cycle (A = 9337, D « 0.57) and the 

other corresponding to N= 1.95* 10 -6 being unstable (A U = 1436, 0 « 1.06). 
These results were not confirmed by analogue simulation. It is felt by the 
author that they are sufficiently straightforward and do not warrant the effort. 


ANALYSIS OF BEND IX MODEL NUMBER 2 


An analysis of Bendix model number 2 was begun. Before it was 
completed, the momentum of opinion appeared to swing strongly in the direction 
of the "Dahl model" (see next analysis section) . The author was swept away 
with this same enthusiasm and terminated analysis of Bendix Model No. 2. 


ANALYSIS OF DAHL MODEL 

The model chosen for the CMG and vehicle single-axis dynamics is 

shown in Figure 4. The nonlinear relationship (indicated by the symbol N) 

between gimbal friction torque T and gimbal angle 6 is portrayed in Figure 

r 

5 and by equations (l) and (2) . As in the above analyses, this analysis was 
based on using a describing function to represent the nonlinearity N . This in 
turn was based on the assumption of a sinusoidal input to the nonlinearity of 
amplitude A and frequency fi. 

Kuo described this nonlinearity by developing a describing function 
[ 11, 12] . Instead of using the straight line approximations conventionally used 
in the derivation of many describing functions, he used the actual curved lines 
precisely representing equations ( l) and (2). Using the complex notation of 
N, 


N = + iN 2 , (49) 
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Kuo obtained the following expressions: 


N 


1 



( 50 ) 


and 


N 


2 



(51) 


where 



4T 


GFo 



(C 2 - A) (C l + A) 
(C 2 + A) (cj - A) 


( 52 ) 


and 


C 2 S 

1 ' ‘ 


B = - — - 

^ yA k f a 2^ * . a kT rZ A? 


( 53 ) 


yA \/ C 9 - A^ yA n/ C * - A^ 


where 


C. = A 


T ^ T GFi " T GFo^ 


( 54 ) 


and 


C 9 = - A - 


y( T~ ‘ T ) 
7V GFi GFo' 


( 55 ) 
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where 


. = _ T 
GFi GFo 


(' a + J 2 ~? i ) 


(56) 


and 


T =■ T 
GFi GFo 




(57) 


where 


a 


2yA T 


GFo 


(58) 


The author found that Kuo’s equations can be simplified considerably by 
observing that C g = - Cj and = -T^., leading to 



■) 

where C = . Alternately these expressions may be rewritten as functions 

of the arguments a and T _ , leading to 

GFo 
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(61) 


A. = 


ln(a+*/^T"i )2 _ 2] 


and 


4T 


B = 
1 


GFo 


a 2 + a+l + (a + l)\/ a 2 + 1 
2(a W a 2 + 1) 


- 1 


(62) 


Approximations for these curves may be obtained by deriving the asymptotes 

for N versus N . The asymptote associated with small values of a (a « l) 
2 1 
is 


N « 7 T 2 
1 7 GFo 


(62) 


and the asymptote associated with large values of a (a » l) is 


N « 
2 


( 4 T GFo ) 


2 / 3 




(64) 


The latter expression is simplified if one plots the real (N^) and imaginary 

(N ) parts of the describing function on log -log paper: 

2 


log N 2 


log K 


+ 3 l ° gN l 


(65) 


where 


K 


^ 4 T GFo^ 


V 3 



(6(5) 
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-«r < 


Hence one obtains a stiaight line with a slope of 2/3 and a y-intercept of K 
for one asymptote and a straight line, parallel to the v-axis, for the other 
asymptote. The effect of variations in the gimbal bearing parameters T 
and y may now be seen explicitly. ° 


Let k T be the factor by which T is increased (or decreased), 
i.e. the new value of T__, is k_ times as large as the former value. Also 

Ur O l 

let k^ be the factor by which y is increased (or decreased). Then, examina- 
tion of equations (63) and (64) shows the effect of altering the numerical 

values of T and y is to multiply the former values of N and N (of 
Or o 12 

the describing functions) by the quantity, kLk . 

T y 


Now, a general map of N 9 versus may be developed for use in 

this and future limit cycle investigations for this particular form of nonlinearity. 

It is universal in the sense that it permits one to choose any value of T 

J GFo 

and y and see its effect on the N versus N map. such as by placing it on 

u 1 

an Nj, N 2 parameter plane stability map to predict and analyze limit cycle 

existence and behavior. 


A new parameter is defined: 



(67) 


Now, a map of N 9 /cr versus N^/o is simply obtained oy plotting the two 
straight lines (asymptotes) of equations (63) ~nd (64) on Figure 32. 

Values of A may also be found from Figure 32 by again resorting to 
approximations. For the asymptote associated with relatively large values of 
a. 


4T 


N * 
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GFo 


7T A 


( 68 ) 



which leads to an expression for \ 
T 

GFo* 


in terms of N <j 


and parameters 


) 


and 


A 


^ T C.Fo^ N 2 


((if) 


2 ' USl " B P<IU;,,i " n -n S c of numerical v^' «'i m 'L 

1.M analysis, equation (no) „ as proven to be a fairly b «hI approximation. 

Let the natural logarithm term (■< ^ i\" » 

K *' n ’ ' ** ' »K be approximated by the term. 


( 


-a 2 iia -* 1 



Then a power series expansion of In f is 



2a(a + l) 
a* T771 + 


2 

3 




a( a + 1) 
a 2 + a + 1 


( 71 ) 


Combining equation (71) with equations (61) and (51) leads to an approximation 
for A as a function of N Jo : 



(72) 


For LST numerical values, this turns out to be a poor, but barely acceptable, 
approximation. The problem in approximations arises, for LST numerical 
parameters, because a difference between two nearly identical numbers is 
required in equation (61) . 

Finally, equations (50), (51), (6l), and (62) are used to determine 
the transition curve connecting the two asymptotes. The describing function 
plot of Figure 32 may now be used in conjunction with a stability contour, such 
as will be developed in the next section, on an N jt parameter plane. If it 

is redrawn on a transparency it may be used with a stability map (as will be 
demonstrated on Figures 34 and 35) without redrawing it each time by dis- 
placing it in both the N and N direction by an amount equal to k? k ; . 

(It has been redrawn on Figures 34 and 35 for the sake of clarity in this paper.) 

As in the previous analyses, the possibility of limit cycle existence will 
be determined. When it is predicted that one exists, its characteristics will 
be determined. 

The model shown in Figure 4 may be described in conventional control 
system form (closed loop transfer function) as 
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(73) 


The characteristic equation associated with this model is 


QK S * 4 ( K „ 4 ir) s> 4 (K i 4 V 3 ’ 4 K . K i s 4 VS ' 0 


(74) 


where 


I1A 

( 

“o I 


(75) 


and 
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1 
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v 


(7ti) 


The system described by Figure 4 appears suitable for describing 
function analysis because it is low pass and the system parameters are assumed 
time invariant. A describing function N is used to represent the nonlinearity 
(herein assumed to be CMG gimbal friction). Recall that for a limit cycle to 
occur, it is required that all characteristic equation (74) roots have negative 
real parts except for one pair which must be purely imaginary roots. This 
condition is determined mathematically and, for ease of visualization, is 
portrayed graphically. In the latter case, two adjustable parameters are 
selected, and N 0 . A correlation between these parameters and the roots 

of the characteristic equation is determined by mapping the stability contours 
from the complex s-plane onto the selected (N^, N o ) parameter plane. For 

the system under consideration, only one such stability boundary exists ami is 
easily found. 
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Analysis of equation (74) indicates the absence on the N , N param- 

1 z 

eter plane of either a stability boundary associated with the real roots of the 
characteristic equation or the boundary at infinity. 

The simultaneous solution of the nonlinear relation of equations (50) 
and (5l) and the purely imaginary root boundary yields the condition for a 
limit cycle, assuming the indicated solution occurs adjacent to a stable region 
(true in this case) . This condition is readily apparent on the parameter plane 
as the intersection (if one occurs) between the purely imaginary root boundary 
and the nonlinear locus of equations (50) and (51) . From this point of inter- 
section the frequency (S2) and magnitude (A) of the indicated limit cycle(s) 
may be determined 3s functions of the characteristics of the nonlinearity and 
of the adjustable parameters. Further, the behavior of each limit cycle when 
a small perturbation is applied to its amplitude and, hence, the nature of 
limit cycle stability also is apparent on the parameter plane. 

Turning to the characteristic equation (74) and setting s = ifi, one may 
obtain the real and imaginary parts of the equation: 


-NjO 2 = - (K 0 Kj - KjS* + I GE fi<) (77) 

and 

-N 2 0 = - (K j Kj - Ktf) , (78) 

which yields the Jacobian 

J = fi 3 > 0 > 0 . (79) 


Solving equations (77) and (78) for the purely imaginary root boundary, 
one obtains 


N 
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i ge b4 - k i !? + k o k i 

j? 


(80) 
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and 


N 


2 



n 2 - 


n 


(81) 


This boundary (cross-hatched curve) is sketched on Figure 33, the nature of 
the Jacobian indicating where the stable region lies. If the nonlinear locus of 
the describing function defined in equations (50), (51), (61), and (62) is also 
drawn on the same N - N 9 plane, it is seen that several conditions can occur, 

depending on the numerical values selected for the system parameters. For a 
fixed set of values for K^, K^, K^, either no intersections (or predicted 

limit cycles) occur (Fig. 34a) or two intersections occur (Fig. 34b). In the 
latter case one limit cycle is stable (with indicated amplitude and frequency 

A S and fi S , respectively) and one is unstable (A U , £? u ) . In that case, if the 
amplitude (A) of the assumed simusoidal input to the nonlinearity is always 

less than A , the output of the nonlinearity will asymptotically approach zero. 

However if A ever exceeds the value A U (such as by an initial condition), 

then A will approach A S (and will approach n S ) . In the limiting case 
where the curves osculate but do not intersect, a single semistable (orbitally 
semistable) limit cycle is indicated. To interpret these phenomena physically, 
one may refer to Figure 4. The input to the nonlinearity is 5. T or describing 
function analysis, it is assumed to be sinusoidal and of the form of equation 

( 3 ) , i« e. , 


5 = A sin fit 


( 82 ) 


The output of the nonlinearity, T^, is then assumed to be of the form 



N5 



( 83 ) 
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■ V + If 6 


= N A sinSH + N A cos fit 

-*■ Ca 


= A 4 + N 2 sin (fit + 0) , <p = tan' 1 (N /N ) . (83) 

(Concluded) 

This is in the form of equation (4) . Using equations (82) and (83) with the 
relations shown in Figure 4, one may now obtain tjje values of the other vari- 
ables (e.g. , 9 and 9) when = 0 and a limit cycle is predicted. 

Typical numerical values of LST were altered slightly, because of the 
passage of time between analyses, and are shown in Table 3 (slightly different 
from Table l). They are described in Reference 13. 

TABLE 3. TYPICAL NUMERICAL VALUES FOR LST 


Parameter 


Numerical Value 


2 x 10 4 ( rad-s) " 1 


3 x 10 3 (rad) -1 


271 Nms (200 ft-lb-s) 
5.012 Nms 2 (3.7 ft-lb-s 2 ) 


1.354 x io 5 Nms 2 (10 5 ft-lb-s 2 ) 


1.354 x io 4 Nm (10 4 ft-lb) 


379 Nms (280 ft-lb-s) 


C.271 Nm (0.2 ft-lb) 


1.477 x ip 5 (Nm rad)* 1 [ 2 x io 5 (ft-lb rad) -1 ) 
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The curve for the purely imaginary root boundary is found by substitut- 
ing numerical values from Table 3 into equations (80) and (81). If Table 3 
values for y and are used for the nonlinear locus representing the 

friction [ equations ( 50) , (5l), (6 1) , and (62)1, it is seen that no intersections 
between the two curves will occur (Fig. 35) . (This is the case portrayed in 
the sketch of Figure 34a.) Hence no limit cycles are predicted and stable 
operation is predicted. 



Figure 35. N - N parameter plane — Dahl model, y = 1.477 x io 5 . 

If the asumptotes are also plotted on Figure 35, one sees that, if during 

the design evolution y and/or T are increased, the case shown in the 

Gr O 

sketch of Figure 34b will occur. In that case two limit cycles are predicted, 
one stable and one unstable. One such example is shown in Figure 36, where 
y is increased to a value of 1.477 x lo fi . The unstable limit cycle has a 
predicted amplitude of 1.0 x io -6 rad and frequency of oscillations of 2.8 rad/s. 
The stable limit cycle amplitude is 2. 1 x io -R rad with a frequency of oscilla- 
tion of 5.2 rad/s. Practically, this means that if A has an initial value that 
is greater than 1.0 x io - * rad it will always converge to a value of 2. 1 x io -R 
(the stable rmit cycle). However, if the value of A is never permitted to 
exceed 1.0 x io _fi rad (the unstable limit cycle), no limit cycle operation will 
occur and the system will be asymptotically stable. 
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>aM* 



Figure 36. N - N parameter plane _ Dahl model 




The above predictions were confirmed with a simulation on an EAI 
231R-V analogue computer, again through the efforts of Mr. P. H. Fisher. 
Figure 37 shows the stable limit cycle results, yielding measured values of 
1.5 x 10 -6 rad and 5.2 rad/s for the amplitude and frequency of oscillation of 
the gimbal angle (6). Figure 38 shows the unstable limit cycle results, with 
values of 2 x 10 ’ 7 rad and 2.4 rad/s for gimbal angle amplitude and frequency, 
respectively. These values are sufficiently close to confirm the predicted 
value within the limits of describing functions analysis assumptions. A sum- 
mary of the analytical predictions and simulation results is given in Table 4. 
Hence for the single nonlinearity considered, i.e. , the output gimbal friction 
relation developed by Dahl, no limit cycle is predicted for presently estimated 
LST numerical parameters. However, if the design characteristics of the LST 

CMGs should allow numerical values of y and/or T substantially larger 

GF o 

than those indicated in Table 3, two limit cycles will occur, one stable, one 
unstable. A more detailed simulation effort has been described by G. S. 

Nurre. 5 

This analysis has considered neither the effect of multiple nonlinearities 
nor the effect of sampling (such as will be performed by an onboard digital 
computer) . A sampled data analysis is now underway to consider the latter. 


CONCLUSIONS 


For numerical values considered to be representative of the LST and its 
CMGs, analysis indicates (and analogue simulation confirms) the absence of 
limit cycle bahavior due to the CMG output gimbal friction nonlinearity. 

Because of the present early stage of development of the LST, it is expected 
that the numerical values (and indeed mathematical characterization of CMG 
friction models) will change. Because of this, a wide spread of numerical 
values of actual variations in friction parameters has been examined. Further, 
a general technique of analysis of such system has been specified and used in 
great detail. This technique should be applicable to future alterations in the 
LST. As such, it provides a design tool for enhancing the efficiency of large 
scale computer simulation by predicting dynamic responses (thus aiding in the 
time-consuming debugging process) and in helping to select numerical values 
to be incorporated in the simulations. While the Dahl model has been widely 
accepted by those involved in CMG dynamics, the author respectfully points 
out that not too long ago the model termed herein as Bendix model number 1 
was widely accepted (at least at MSFC) . 


5. G. S. Nurre, An Analysis of the Dahl Friction Model and its Effect on a 
CMG Gimbal Rate Controller, S&E^ASTR-A, Unnumbered report, MSFC, 
Oct. 18, 1973. 
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Figure 37. Dahl model simulation, y * 1.477 x 10' . 


'l'he analysis of Hendix model number 1, in particular, is quite detailed, 
perhaps to the point of inciting boredom. However, not only is the general 
analysis technique presented but also a study of the effect of parameter varia- 
tions is presented. Also, from an academic viewpoint, it is interesting (to 
some) to see actual examples of such usually obscure limit cycle behavior as 
orbitally semistable limit cycles. 
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